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Abstract 

We consider the problem of computing the minimum of a polynomial function g on a 
basic closed semialgebraic set E C M". We present a probabilistic symbolic algorithm to find 
a finite set of sample points of the subset E^™™ of E where the minimum of g is attained, 
provided that E™^'^ is non-empty and has at least one compact connected component. 
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"S ' 1 Introduction 

S ■ 

The ininimization of polynomial functions over R", unrestricted or subject to polynomial con- 
straints, is a classical problem with a variety of applications. In the last years, it has been exten- 
sively studied in the algorithmic framework through numerical or symbolic-numerical methods 
^ ; based on certificates of positivity (see, for instance, [20], [22], [23], [21], [27], [II], |9]). 

l/-^ ■ From the symbolic computation perspective, a possible way to tackle the problem is to restate 

•^ , it as a quantifier elimination problem over the reals and to apply a symbolic algorithm to solve 

this more general task. However, even the most efficient algorithms for quantifier elimination (see 
^T ' [U Chapter 14]) lead to high complexities since they are not particulary designed for this specific 

^f^ . problem. In [25], a probabilistic algorithm for global optimization of polynomial functions over 

M" with good complexity estimates is presented. The better complexity is due to an alternative 
strategy relying on the computation of generalized critical values. Also, the problem of deciding 
algorithmically whether the global infimum of a polynomial function is attained is considered in 

m- 

Recently, in [H] (see also [6] , [H] ) , another approach based on computer algebra techniques 
led to a lower bound for the minimum of a polynomial function on a basic closed semialgebraic 
set, provided that the set where this minimum is attained is compact. Here, we adapt the 
techniques underlying these theoretical result to the algorithmic framework. 

More precisely, we consider the following problem. Let K C M be an effective field. Let 
fi,---,fm,g^ K[xi,...,Xn], with n > 2, and 

E = {x E M" I /i(x) = • • • = fi{x) = 0, fi+i{x) > 0, . . . , U{^) > 0}, 
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and suppose that g attains a minimum value 5min at E. We look for a symbolic algorithm to 
compute at least one point in 

E^'^ = {xeE\ g{x) = <7min}. 

In this paper we assume that E"^™ has at least one compact connected component. This is 
the case, for instance, in many families of instances of known optimization problems (see |18t 
Section 4]). Our approach consists in finding a finite set of points containing at least one point in 
each compact connected component of i?™"'. The natural tool to use when solving this problem 
is the Lagrange Multiplier's Theorem; nevertheless, a direct application of this result may lead 
to a degenerate system or a system with infinitely many solutions. In order to overcome this 
difficulty, we apply deformation techniques as in [T7] (see also [H Chapter 13]), which enable us to 
deal with "nice" systems that, in the limit, define finite sets containing the required minimizing 
points. These sets are described by geometric resolutions, which are parametric representations 
where the parameter ranges over the set of roots of a univariate polynomial. Finally, we compare 
the values that the given polynomial function g takes at the computed points and obtain the 
Thom encodings characterizing the minimizers. 

The main result of the paper is the following: 

Theorem 1 Let E = {x eW \ fi{x) = ■ ■ ■ = fi{x) = 0, fi+i{x) > 0, . . . , /^(a;) > 0} he defined 
by polynomials /i, . . . , /m £ ^[xi, . . . , x„] with n > 2 and degrees bounded by an even integer 
d. Let g E K[xi, . . . ,3;„] be a polynomial of degree at most d that attains a minimum value 
5min at E in a non-empty set E)™™ with at least one compact connected component. Algorithm 
FindingMinimum (see Section \^\j is a probabilistic procedure that, taking as input the integer d 
and the polynomials fi, ■ ■ ■ , fm,g encoded by a straight-line program of length L, computes a 
family 

where X is a finite set and for every i S Z, {pi,Vi^i, . . . ,Vi^n) is a geometric resolution in K[ii] 
and Ti S { — 1,0, Ij'^^SP' is the Thom encoding of a real root ^j of pi such that the set 



{{viAii), ■ ■ ■ ,Vi^n{ii))}; 



GX 



is included in i^;™™ and intersects all its compact connected components. The complexity of the 
algorithm is 



0[ {n\L + dn + n^^~^)D'^ log%D) log log^(i:)) + (m + D)D'^ log\D))T 



where 




• D = max (' 

0<s<min{n,m} 


:)d'id-l 


0<s<niin{n,m} 


E 

SI + S2 — S 



0<s<min{n,m} 

In the above statement, O denotes a positive real number such that for any ring R, addition, 
multiplication and the computation of determinant and adjoint of matrices in R''^'' can be 
performed within 0{k ) operations in R. We may assume Q < 4 (see [3]) and, in order to 
simplify complexity estimations, we also assume that fi > 3. 



An extended abstract of this work containing Theorem [T] has been accepted for presentation 
at the conference Effective Methods in Algebraic Geometry (MEGA) 2013 (|16]). 

The paper is organized as follows. In Section [21 we introduce the basic notation and state 
some previous results we use throughout the paper. In Section [31 we present the deformation 
we apply and we prove some of its geometric properties. Section [D is devoted to showing how 
the deformation leads to a geometric resolution of the finite set we look for. Finally, in Section 
[5] we show the algorithmic counterparts of the previous theoretical results, proving Theorem [H 

2 Notation and preliminaries 

Throughout the paper, we denote N the set of positive integers. For a field K, we write K for 
an algebraic closure of K^ K{t) for the field of rational functions in a single variable t and i^[[t]] 
for the set of formal power series in t. 

For a given n € N, the n-dimensional affine and projective spaces over an algebraically closed 
field K are denoted by A^ and P^ respectively. In the case when i^ = C, we write simply A" 
and P*". 

2.1 Algorithms and complexity 

The algorithms we consider in this paper are described over an effective field K C M. The notion 
of complexity of an algorithm we consider is the number of operations and comparisons over K 
that the execution of the algorithm requires. In this definition of complexity, accessing, reading 
and writing pre-computed objects is cost free. 

The objects we deal with are polynomials with coefficients in K. In our algorithms we 
represent each polynomial either as the array of all its coefficients in a pre-fixed order of its 
monomials {dense form) or by a straight-line program. Roughly speaking, a straight-line program 
(or sip, for short) over K encoding a list of polynomials in K[xi, . . . ,Xn] is a program without 
branches (an arithmetic circuit) which enables us to evaluate these polynomials at any given 
point in K". The number of instructions in the program is called the length of the sip (for 
a precise definition we refer to [H Definition 4.2]; see also [l3]). From the dense encoding of 
a family of m polynomials in K[a;i, . . . ,Xn] of degrees bounded by d, we can obtain an sip of 
length 0(771 ( ^")) encoding them. Unless otherwise stated, throughout the paper, univariate 
polynomials will be encoded in dense form. 

To estimate complexities we will use the following results. Operations between univariate 
polynomials with coefficients in a field K of degree bounded by d in dense form can be done 
using 0((ilog((i) loglog((i)) operations in K (see [71 Chapters 8 and 9]) and gcd or resultant 
computations by means of the Extended Euclidean Algorithm (see [71 Chapter 11]) can be 
performed within 0(dlog {d)loglog{d)) operations in K. Given an sip of length L encoding 
a univariate polynomial of degree at most d, we can obtain its dense form within 0{dL + 
(ilog^((i) loglog((i)) operations in K (see [3 Corollary 10.12]). 

From an sip of length L encoding a polynomial / € K[xi, . . . , Xn], we can compute an sip of 
length 0{L) encoding / and all its first order partial derivatives (see [2]). 



2.2 Geometric resolutions 

A way of representing zero-dimensional affine varieties which is widely used in computer algebra 
nowadays is a geometric resolution (see, for instance, [5]). The precise definition we are going 
to use is the following: 

Let K he & field of characteristic and V = {zi, . . . ,zd} C A^ be a zero-dimensional variety 
defined by polynomials in K[xi, . . . , Xn]- Given a separating linear form i = aixi + • • • + anX„, G 
K[xi, . . . ,Xn] for V (that is, a linear form I such that i{zi) ^ ^(zj) if i ^ j), the following 
polynomials completely characterize the variety V: 

• the minimal polynomial p := ni<i<D(^ ~ ^(-^0) ^ K[u] of i over the variety V (where u 
is a new variable), 

• polynomials vi, . . . ,Vn G K[u] with deg{vj) < D for every 1 < j < n satisfying 

v = {{vi{0,..., ^n(O) (^K"\^eK, piO = 0} . 

The family of univariate polynomials {p,vi,. . . ,f„,) is called a geometric resolution of V (asso- 
ciated with the linear form i). Note that if X is a subfield of M, the real roots of p correspond 
to the real points of the variety V. 

Given geometric resolutions (pi,fii, . . . ,vin) and {p2,V2i, ■ ■ ■ ,V2n) of two zero-dimensional 
varieties Vi and V2 in A" consisting of Di and D2 points respectively, associated with the same 
linear form i which separates the points in Vi U V2, we can obtain a geometric resolution of 
Vi U V2 within complexity 0{nDlog (D) loglog(L')), where D = m.ayi{Di, D2}, by means of the 
Chinese Remainder Theorem using the Extended Euclidean Algorithm. 

2.3 Thorn encoding of real algebraic numbers 

The Thorn encoding of real algebraic numbers provides an algebraic approach to distinguish 
the different real roots of a real univariate polynomial. We recall here its definition and main 
properties (see [TJ Chapter 2]). 

Given p € K[u] and a real root ^ of p, the Thom encoding of ^ as a root of p is the 
sequence (sign(p'(^)), . . . ,sign.{p^'^^^P'{^))), where we represent the sign with an element of the 
set { — 1,0, 1}. If the sign of the leading coefficient of p is known, the Thom encoding can be 
shortened to (sign(p'(^)), . . . , sign{p^'^'^^P~^\^))) . 

Two different real roots of p have different Thom encodings. In addition, given the Thom 
encodings ((Ti^i, . . . , (Ti^dcgp) and (o"2,i, . . . , '72,dcgp) of two different real roots ^1 and ^2 of p, it is 
possible to decide which is the smallest between ^1 and ^2 as follows: Consider the largest value 
of k such that ai^k ¥" '^2,k', then k < degp, since p('^°sp) jg g, constant. Also ai^k+i = (^2,k+i 7^ 0, 
since otherwise ^1 and ^2 would have the same Thom encoding with respect to the polynomial 
p(«+i) and therefore ^1 = ^2- Then, 

• if CFi^k+i = o-2,fc+i = 1; we have that ^1 < ^2 if and only if ai^u < cr2,k, 

• if (Ti^fc+i = cr2^k+i = —1) we have that ^1 < ,^2 if and only if ai^k > (^2,k- 



3 The deformation 

3.1 Defining the deformation 

Here we introduce the deformation we use. We denote: 

• qo := n + 1 and qi < ■ ■ ■ < Qm the first m prime numbers greater than n + 1. Let 
A G qH+i)x ("■+!) J ^ = {aij)o<i<rn, o<j<n be the Cauchy matrix defined by aij = —^ 
(note that each submatrix of A has maximal rank and a^j > for every i,j). 

• For e S N, Te the Tchebychev polynomial of degree e (see |19t Section 6.1]). 

• 9{^) = Y^i<j<n^OjTd{xj) and, for every I < i < m, fi{x) = ajo + Y.i<j<n^ij{Td{xj) + l). 

• G{t, x) = tg{x) + (1 — t)g{x) and, for every 1 <i <m^ ^ti^^ ^) ~ ^/K^) + (1 ~ t)fi{x) ^iid 
Fr{t,x)=tfi{x)-{l-t)h{x). 

• For every 5 C {1, . . . , m} and a G {+, — }'^, 

Vs,a = { (t, X, A) G A X A" X Pl'^l I Fp {t, x) = for every i G S, 
AoV,.G(t,x) = J^AiV,.i^'^'(t,:E) }. 

We consider the decomposition of Vs,a as Vs,a = ^si^ ^s,cr, where 

— y^ ^ is the union of the irreducible components of Vs^a included in i = io for some 

to e C, 

— Vs,a is the union of the remaining irreducible components of Vs,a- 

• For a group of variables y, Hy the projection to the coordinates y. 

3.2 Geometric properties 

Let C be a compact connected component of E"^"^. For 5 > 0, we write 

• C=s = {x G R" I dist(x, C) = 6}, 

• C<s = {x G M" I dist(x,C) < 6}, 

• C7<5 = {x G M" I dist(x, C) < 6}. 

Let /i > such that C<^ and -Emin \ C do not intersect. 
We consider 

E = {{t,x) gMxR" I F+(i,x) >0, ...,F+(t,x) > 0, F[^^{t,x) > 0, . . . , F+(t, x) > 0, 

F^{t,x)<0,...,Fr{t,x)<0} 

and, for every < t < 1, 

Ft = U^{En{{t} X M")) C M". 

Note that Fq = M", Fi = F and for < ti < ^2 < 1, Ft^ C Ft^^. 



Lemma 2 There exists < e < 1 such that for every 1 — £ < t < 1, the minimum value that 
G{t, ■) takes on Et n C<^ is not attained at any point in Et n C=^. 

Proof. The minimum value oi g{) = G{1,-) at Eif] C<^ is gmim which over this set is only 
attained at points in C and, therefore, not in Ei r\C=^. Assume the claim does not hold. Then 
there exists a strictly increasing sequence of positive numbers {tk)k&n converging to 1 such that 
the minimum value of G{tj^, •) at Et,, n C<^ is attained at a point z^ € Et^. n C=^. 

The sequence {zk)k&i is contained in the compact set C=^; therefore, without loss of gener- 
ality, we may assume that it converges to a point z G C=^. On the other hand, the sequence 
(ife)'2fc)fc£N is contained in E and converges to (l,^;), then (1,2) € E and z G Ei. 

We will prove that g{z) = ^min- Take z' € Ei n C<^, then z' G Et,^ fl C<^ for every A; G N 
and 

g{z') = G(l, z') = hm G{tk, z') > lim G{tk, z^) = G(l, z) = g{z). 

Then g attains its minimum on Ei nC<^ at z, which is grain- This leads to a contradiction since 
this value is not attained at any point in Ei n C=^. D 

The following proposition shows that in order to obtain minimizers for the polynomial func- 
tion g on the compact connected component C, it is enough to consider at most as many of the 
equations and inequations defining E as the number of variables. We define the set 

S = {(5,cr) I S C {!,... ,m} with < |5| < n and a G {+,-}"^ with Oi = + ioi l + l < i < m] . 

Proposition 3 Let C he a compact connected component of E"^^^ . There exist z ^ G and 
{S, a) eS such that z G 112,(^5,^ D {t = 1}). 

We point out that if E (resp. E'™™) is compact, we can easily prove Proposition [3] adapting 
the arguments in the proof of \lii\ Proposition 7] (resp. |18| Theorem 14]). Now we prove it 
under our weaker assumptions. 

Proof of Proposition [21 Consider < e < 1 such that 

• for every 1 — £ < t < 1, the minimum value that G{t, •) takes on Et fl C<^ is not attained 
at any point in Et n C=^, 

• for every S" C {1, . . . , m} and cr G {+, -}^, nt(yj*]) n (1 - e, 1) = 0, 

• for every 1 — e<i<l, S'C {1, . . . , m} with l^l > n and a G {-|-, —}^, the set 

{x G A" I F[' [t, x)=0 for every i G S} 
is empty. 

The existence of such an e follows from Lemma [21 the finitness of '^t{Vs i) ^°^ every S and 
o", and an adaptation of the arguments in \n\ Lemma 21] or [181 Lemma 4]. 

Let (tfc)fcGN be an increasing sequence converging to 1 with ii > 1 — e and let Zk G Et^ fl C<^ 
be a point such that G{tk, •) attains its minimum value on the set Et,. nC<^ at Zk- Without loss 
of generality, we may assume that the sequence {zk)keN is convergent to a point z G Ei n C<^, 
and proceeding as in the proof of Lemma [21 we have that z (z C. 



Now, for every A; € N and every x € M", at most one F^{tk,x) and F^ {tk,x) may vanish. 
Let 

Sk = {ie{l,...,l}\ F+{tk, Zk) = or Fr{tk, Zk) = 0} U {i e {I + I, . . . ,m} \ F+itk, Zk) = 0}; 

then < \Sk\ < n. Without loss of generahty, we may assume that Sk is the same set S for 
every A; € N; moreover, we may assume that, for each i G S (1 {1, . . . ,1}, it is always the same, 
F^{tk,Zk) or F,[~{tk, Zk), the one which vanishes, thus defining a function a € {+,— } with 
(Tj = + for / + 1 < i < m. Then, (5, a) £ S. 

For A: € N, if the set {VxF^^ {tk, Zk),i G S} is linearly independent, since the function G(tfc, •) 
attains a local minimum at the point Zk when restricted to the set Et^ n C<^ , by the Lagrange 
Multiplier's Theorem, there exists {\i^k)i<^s such that 

VxG{tk,Zk) = y2Xi,k'^xFp{tk,Zk). 

We take Ao,fc = 1 and we have that {tk,Zk, (Ao,fc, {Xi,k)ies)) S Vsy, but since tk nj(y^*]), we 
have that {tk, Zk, {Xo,k, {Xi,k)ies)) G Vs^a- On the other hand, if Y^i^^s Xi,k'^xF['{tk,Zk) = with 
(Ai,fc)ies / 0, we take Ao,fc = and, as in the previous case, (tk, Zk, (Ao,fc, {Xi,k)i€s)) S ^5,^^- 

Without loss of generality, we may assume that (Ao,fc, (Ai,fc)ies)feGN converges to a point 
(Ao, {Xi,o)ies)) e P'*^'; then {l,z, (Aq, (Ai,o)je5)) e Vs,a and, therefore, z € Ux{Vs,a n {t = 1}) as 
we wanted to prove. D 

Proposition 4 For every {S,a) € S, we have that Ilx{ys,<j H {t = 1}) is a finite set. 

The proof of Proposition [3] will follow from arguments in the next section. From Propositions 
E] and m we deduce the following: 

Corollary 5 Under the assumptions of Theorem\^ the set 

U l[x{Vs,.r\{t = i}) 

is,a)es 

is finite and contains a point in every compact connected component of E™^'^. 

Due to Corollary [5l to solve the problem we are considering, in the next section we will focus 
on describing the set Tixiys,a n {t = 1}) for every {S,a) G S. 

4 A geometric resolution 

Throughout this section, we consider fixed (-S, cr) e S and denote V = Vs^a and V = Vs,cr] 
moreover, for simplicity, we suppose S = {1, . . . , s} and cr = {+}^ . 
For 1 < j < n, let 

9j{x,X) = Xq— - ^ Aj— ^ G K[xi,...,x„,Ao,...,As], 

-' l<j<s ■' 



gj{x,X) = Xq— - ^ Aj— -^ G K[xi,...,x„,Ao,...,As], 

Gj{t,x,X) =tgj{x,X) + (1 - t)^j(x,A) G K[t,xi,... ,Xn, Aq, . . . , A^]. 

These polynomials are homogeneous of degree 1 in the variables A; therefore, by the multi- 
homogeneous Bezout theorem (see, for instance, |26t Chapter 4, Section 2.1]), the degree of the 
varieties V D {t = to} for to S C is bounded by 

The next lemma shows the key properties of the initial system in the deformation (c.f. [171 
Lemma 20]). 

Lemma 6 The polynomials fi, . . . , fs,gi, ■ ■ ■ ,gn define a 0- dimensional variety in A" x {Aq 7^ 
0} C A" X P* with Ds distinct points wi, . . . , wd^ satisfying Yix{wi) 7^ Jlx{wj) for i 7^ j, and the 
Jacobian determinant of fi, . . . , fs and the polynomials obtained from gi, . . . ,gn dehomogenizing 
with Ao = 1 does not vanish at any of these points. 

Proof. Recalling that g{x) = I]i<j<nOoj7d(xj) and fi = Ojo + Z]i<j<n"ii(^rf(^j) + 1)' '^e have 
that, for every 1 < j < n. 



gj{x,X) = T'^{xj){aojXQ- ^ ajjA, 



l<i<s 

Therefore, the solution set of the system fi, . . . , fs,gi, . . . ,gn can be decomposed as 
U {r^(x,) = Vj G B, h{x) = 0, . . . , fs{x) = 0} X {ao.Ao - Y. «^i^^ = V j ^ B}. 

-Bc{l,...,n} l<i<s 

By our assumption on the matrix A,\{\B\=n — s, the linear system oqjXq — X^i<j<g a-ijXi = 
yj ^ B has a unique solution A^ G P''; moreover, this solution lies in {Aq 7^ 0}. 

For a fixed i? C {1, . . . , n} with \B\ = n — s, taking into account that T^ has d — 1 real roots 
and Td takes the value 1 or —1 at each of these roots, we have that 

Sb := {r^(x,) = Vi G B, h{x) = 0, . . . , fs{x) = 0} 

decomposes as the union of the sets 

SB,e ■■= {T'Axj) = 0, Td{xj) = e{j) Vi G B, /f '^ = 0, . . . , /f '^ = 0} 

for all e : B — ?• {1,-1}, where /^ ' G K[xj; j ^ B] denotes the polynomial obtained from fi by 
replacing T(i{xj) = e{j) for every j £ B. 

Without loss of generality, in order to simplify notation, assume B = {s + 1, . . . ,n}. Then, 
for e : B ^ {1, —1}, the system f^ ''^ = 0, . . . , fs '^ = can be written in the form 

B,( 



B,e 

a. 




Ci 





° ^ 


* 





C2 


* 


C3 


* / 



where Ab := {aij)i<i,j<s and a- '*" = -ajo - I]s+i<j<n "^^iC^O) + 1) for 1 < i < s. Since 
^4^ is invertible, we can solve the underlying linear system for Td{xi) + 1, . . . ,Td{xs) + 1. By 
applying Cramer's rule, it can be seen that the coordinates of the solution to this linear system 
are rational numbers where the denominators are a multiple of the prime number q^, whereas 
the numerators are relatively prime with Qs] therefore, no coordinate of a solution is an integer 
number. We deduce that the above system is equivalent to a system of the form 

Td(xl) = cf'^...,^,(x,) = cf•^ 

where c^ ' 7^ ±1 for every 1 < i < s. It follows that each of the equations has d distict roots, 
none of which equals a root of T^ (thus, the sets 5^,6 are mutually disjoint). 

Moreover, the Jacobian matrix of fi, ■ ■ ■ , fs and the polynomials obtained from gi, . . . ,gn 
dehomogenizing with Aq = 1 evaluated at any of its solutions is of the form 

s { 
s { 

n — s { 



It is easy to see that Ci, C2 and C3 are invertible matrices and so, the Jacobian determinant 
does not vanish. 

We conclude that Sb consists of {d—l)^~'^d^ distinct points in A" for every B with \B\ = n—s. 
Hence, the system fi, . . . , fs,gi, . . . ,gn has Ds isolated solutions in A" x P* whose projections to 
A" are all distinct. Since Dg is an upper bound for the degree of the variety the system defines, 
it follows that these are all its solutions. D 

As a consequence of Lemma El it follows that all the irreducible components of V intersect 
the set {t = 0} and have dimension 1. Proposition H] is immediate from this fact. Moreover, the 
following further properties of the induced deformation hold. 

Lemma 7 The variety defined in A^^^-— - x P^— ry by Fi, . . . ,Fs,Gi, . . . , Gn is 0- dimensional and 

has Ds distinct points Wi, . . . , Wds in {Aq 7^ 0} such that HxiWi) ^ Ilx{Wj) for i ^ j- More- 
over, these points can be considered as elements in K[[t]]"~^'*. 

Proof: The multihomogeneous Bezout Theorem states that the degree of the variety is bounded 
hy Ds- If wi, . . . , wDs are the common zeros of fi, . . . , fs,gi, ■ ■ ■ ,gn, the Jacobian with respect 
to xi, ... ,Xn,Xi, ... ,Xs of Fi, ... ,Fs, and the polynomials obtained from Gi, . . . , G„ dehomog- 
enizing with Aq = 1 at t = and (x. A) = Wi is nonzero. The result follows applying the 
Newton- Hensel lifting (see for example [12^ Lemma 3]). D 

Consider now new variables yi,...,y„ and define £{x,X,y) = i{x,y) = J2i<j<nyj^j- ^^^ 
ai, . . . , Q„ € C, let iaix, A) = i^ix) = I]i<j<n "j^j- Let 

r>u \ TT /' fffui \\ T.o<h<DsPh{t:y)u P{t,u,y) 
P{t,u,y)= [[ [u-i{W^,y)} = - , . = -7 — €K{t)[u,y], 

l<i<D. ^^ ^ ^^ ^ 



with P(t,u,y) € K[t, n, y] with no factors in K[f]. Let Q{u,y) = gcd{P {I, u, y) , dP / du{l, u, y)) . 
Then, for generic a £ C", if 

P{l,u,a) 



p{u) :-- 



Q{u,a) 



and, for every 1 < j < n. 



Q{u,a) \ Q{u,a) I 

we have that (^p{u),vi{u), . . . ,Vn{u)^ is the geometric resolution associated to £a of a finite set 
V containing na;(y {t = 1}) (see, for instance, [SI Algorithm 9] or [171 Proposition 8]). 

The computation of P will be done by means of a Newton-Hensel based approximation. The 
required precision is obtained from the degree bound deg^ P{t, u, y) < nDg, which can be proved 
as in [171 Lemma 9]. 

5 The algorithm 

In this section we present our algorithms and prove the main result of the paper. 
From Corollary [SJ we know that the finite set 

U u,{Vs,^n{t = i}) 

iS,a)£S 

contains a point in every compact connected component of £^™™; nevertheless, for a fixed (S, a) S 
S, the set Ilxiys,<j H {i = 1}) is not necessarily contained in i?™"^, or may even have an empty 
intersection with E. The idea of our main algorithm is to compute first finite sets 1^3,0 containing 
Ilx{ys,a n {t = 1}); then, look for for the points of each 'Ps,a that lie in E and finally, compare 
the values that the function g takes at these points. 

First, we introduce three auxiliary subroutines we use to construct our main procedure. 

Our first subroutine is an algorithm to compute the geometric resolution of the finite set 
Vs,a containing IlxiVs^a n {t = 1}) introduced in the previous section. This algorithm relies on 
the global Newton lifting from [8] and it is essentially the procedure underlying |17^ Proposition 
13]; we include it here for the sake of completeness. In order to simplify notation, we assume 
that 5 = {1, ..., s} and cr = {+}^ ■ 

Algorithm GeometricResolution 

INPUT: Polynomials /i, . . . , /s, g G K[xi, . . . , Xn] encoded by an sip of length L, an even integer 
d > deg{fi),deg{g), and a linear form ia € K[xi, . . . , Xn]- 

OUTPUT: The geometric resolution {p,vi, . . . ,Vn) associated to ia of a finite set V containing 

Ux{vn{t = i}). 

1. Compute the geometric resolution associated to ia{x) = aixi + • • • + anXn of the variety 
defined in A"+* by the (dehomogenized) system fi, . . . , fs,gi, . . . ,gn as follows: 



10 



(a) For every B C {1,... ,n} and e : B —^ {—1,1}, compute the geometric resolution 
associated to £aix) of the variety SB,e- 

(b) Compute the geometric resolution associated to £a{x) of the variety |J^ ^ SB,e- 

2. Compute the geometric resolution associated to i{x, y) = yixi + • • • + ynXn of the variety 
defined by the (dehomogenized) system fi, . . . , fs,gi, ■ ■ ■ ,gn over K(y), modulo the ideal 
(yi -ai,...,y„ - a^)^. 

3. Compute P{t, u, y) mod {{tf^D,+^ + (y^ _ a^, . . . , y„ - a„)2)K[[t]][n, y]. 

4. Compute P{l,u,a) = Y.o<h<D,Ph{'^,oi)u'^ and ^(l,u,a) = Eo<ft<D, ^ (!>«)«'' as 
follows: 

(a) Compute ph{t, a) and ^{t,a) (l < j < n, < h < Ds). 

(b) Evaluate t = 1. 

5. Compute Q{u,a) = gcd(P(l, n, q), ^{l,u,a)), p{u) = ^^^^ and, for every 1 < j < n, 

Proposition 8 Given a generic a G K" anrf polynomials fi, ■ ■ ■ , fm G K[xi, . . . , 3;„] of degree 
bounded by an even integer d and encoded by an sip of length L, Algorithm. GeometricResolution 
computes the geometric resolution associated to the linear form ^q(x) = X^i<,<n'^j^j of a finite 
setVs,a containing Ilx{Vs^a(~^{i = 1}) within complexity 0{n^{L+dn+n^^)Dl\og (D^) log log (Ds))- 

Proof: Step 1(a). The variety SB,e is defined by a square polynomial system in separated vari- 
ables; then, the required computation can be achieved as in |14| Section 5.2.1] within complexity 
0(D|glog^(L's,e)loglog(DB,e)), where Db^b is the cardinality of Sb^c- 

Step 1(b). This step can be done within complexity 0{nDs\og {Ds)\og\og{Ds)) following the 
procedure in Section YI?2\ and the strategy described in [71 Algorithm 10.3]. 

Step 2. This step can be done applying [HI Algorithm 1] within complexity 0{{dn^ + n ^^) 
Ds\og{Ds)\og\og{Ds)). 

Step 3. Since Fi, . . . , Fs, Gi, . . . , Gn can be encoded by an sip of length 0{L + (d + s)n), 
a geometric resolution of the variety they define associated with the linear form i{x, y) mod- 
ulo the ideal (i)2nDs+i _|_ ('y^ _ q,_^^ _ _ _ ^y^ _ q,,^)2 ^^^^ |^g obtained from the previously com- 
puted geometric resolution by applying [HI Algorithm 1] within complexity 0{n^{L + dn + 
n^-^)Dl\og\D,)\og\og\D,)). 

Step 4(a). By expanding P{t,u,y) = '}2a<h<DsqW"^^ ^ K[[t]][u,y] into powers of u, 
{yi — ai), . . . , {yn — Q-n), we have that the coefficients corresponding to u^ and u^{yj — aj) 
(1 < j < n, < /i < Dg) are ph{t,a)/q{t) and ■^{t,a)/q{t) respectively. As the degrees 
of the polynomials involved in these fractions are bounded by nDs, they are uniquely deter- 
mined by their power series expansions modulo (t)2"^=~'"^K[[t]] (see [3 Corollary 5.21]) that 
were computed at Step 3. By using [TJ Corollary 5.24 and Algorithm 11.4] and converting 
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all rational fractions to a common denominator, the computation is done within complexity 
0{n^D^,log^{Ds)loglog{Ds)). 

Step 5. This step is achieved by means of the Extended Euclidean algorithm and polynomial 
divisions with remainder within complexity 0{nDs log (Dg) loglog(Ds))- D 

The second subroutine presented here tells us, for a finite set Vs,a given by a geometric 
resolution, if the set Vs,a H E is empty, and, if not, it computes the Thom encoding of all the 
real roots of the minimal polynomial ps.a corresponding to the points where the minimum value 
of g on Vs.cT ri E is attained. 

Algorithm MinimumlnGeometricResolution 

INPUT: A geometric resolution {ps,a,vs,a,i, ■ ■ ■ ,vs,a,n) in K[u] of a finite set Vs,ct, polynomials 
/i, . . . , fm, 5 € K[xi, . . . , Xn] encoded by an sip of length L, an integer < I < m and an integer 
d> deg{fi),deg{g). 

OUTPUT: A boolean variable "Empty" with the truth value of the statement "The set 'Ps,af^E 
is empty" and a list of elements ti, . . . ,rfc € {—1,0, ly^'^sps,^-^ with A; = if Empty = True, 
representing the Thom encodings of all the real roots of ps,a corresponding to the points where 
the minimum value of g on Vs^a n -E is attained. 

1. Compute the list of realizable sign conditions for fiivs^aiu)), . . . , fm{vs,cT{u)) over the real 
zeros of ps^a{u)- 

2. Determine Empty going through the obtained list of realizable sign conditions. 

3. If Empty = False: 

(a) Compute h{u) := ReSu{ps,aiu),u - g{vs,aiu))). 

(b) Compute the list of realizable sign conditions for fi{vs,cTiu)), ■ ■ ■ , fm{vs,a{u))-, p's ai'^)-' 
...,p''s7^^'^~^\u), h'{g{vs,a{u))), . . . ,h^'^^^Ps,--^\g{vs,a{u))) over the real zeros of 
PsA^)- 

(c) Determine ti, . . . , r^ going through the obtained list of realizable sign conditions. 

Proposition 9 Given a geoemetric resolution {ps,(T,vs,a,i, ■ ■ ■ ,vs,(T,n) of the set Vs,a and the 
polynomials fi, . . . , fm,g G K[3;i, . . . ,x„] encoded by an sip of length L, Algorithm 
MinimumlnGeometricResolution decides whether T's,a T\ E is empty or not and computes the 
Thom encodings of the real roots of ps^a corresponding to the points where the minimum value 
of g on Vs,a Pi E is attained within complexity 0{nDl + LDg + (m + Ds)D^ log (Ds))). 

Proof. Step 1. First we compute the dense encoding of the polynomials fi{vs,aiu)), ■ ■ ■ , fm{vs,a{u)), 
within complexity 0{{nDs + L)Ds+mDs log^(Ds) log log(Ds)). Then we apply the sign determi- 
nation algorithm from [U Section 10.3] and [S], following [23], within complexity 0{mD^ log (Ds))- 

Step 2. This step can be done within complexity 0{mDs). 

Step 3(a). At this step the algorithm computes the monic polynomial h whose roots are the 
values of g at the points in 7^3,(7- First we compute the dense encoding of the polynomial 



12 



g{vs,a) and then, the resultant polynomial by multi-point evaluation and interpolation within 

complexity 0{{nDs + L)D, + Dl \og^{Ds) log log(i?,)). 

Step 3(b). We continue the sign determination algorithm adding the polynomials p's„{u), 

■■■^Psa ^'"^ (^)' ^'i9i''^s,(7{u))), . . . ,h^'^^^P^-'^~^'{g{vs,a{u))) to what we have already com- 
puted at Step 1. To do this, we first obtain the dense encoding of the polynomials involved. 
This step can be done within complexity 0{{m + Ds)D^ log (-D^)). 

Step 3(c). The list of sign conditions computed at Step 3(b) enables us to know the Thom 
encoding of every real root ^ oi ps,(T-, and to relate ^ with the Thom encoding of g{vs,a{0) ^^ ^ 
root of h. This information is enough to compare the different values of g{vs,aiO) (sse Section 
12. 3p and, so, we can give the Thom encoding as roots of ps,a of the roots giving the points where 
the minimum value of g is attained. This step is done within complexity 0{{m + Ds)Ds). □ 

Since we know that the minimum value of g on i? is also the minimum value of g on 

U n,.(y5,<xn{t = i})n^, 

{S,a)(iS 

the set of sample minimizing points will be obtained by comparing the minimum values that g 
takes on Vs,(t H E for all (S", a) G S. This task is achieved by the following subroutine. 

Algorithm ComparingMinimums 

INPUT: Geometric resolutions {psi,ui,vs^m^,i, ■ ■ ■ ,vs^,ai,n) and {ps2,u2^vs2,a2,i^ ■ ■ ■ :Vs2,u2,n) of 
finite sets T'si,ai and Vs2,a2 associated with a linear form i^, g G K.[xi, . . . ,Xn] encoded by 
an sip of length L, an integer d > deg{g), and Thom encodings ri S { — 1,0, l}<^'^gPSi,<Ti-i ^nd 
T2 € {—1,0,1} °gPS2,<T2-i Qf j.gal roots oi psi,ai and PS2,a2 corresponding to points where the 
minimum of g on Vsi,ai H E and Vs2,ct2 H -E is attained. 

OUTPUT: An integer "Sign" from the set { — 1,0, 1} representing the sign of the minimum value 
that g takes on Vsi,ai H E minus the minimum value that g takes on Vs2,(T2 l^ -^• 

1 . Compute a geometric resolution {p,vi, . . . ,Vn) of the union of the sets described by {ps^ ,0-1 , 

VSi,ai,l,- ■ ■ ,VSi,cn,n) and {pS2,a2,VS2,a2,l, ■ ■ ■ ,VS2,a2,n)- 

2. Compute h{u) := KeSu{p{u),u — g{v{u))). 

3. Compute the list of realizable sign conditions for p5j^CTi('w) 5 p'sj 0-1 (^)' ••• '^5 o- ^ "^^ (^)' 
PS2,a2(w),Ps'2,<72W' • • • 'Ps2*'^r'''"' N' f^'isiviu))), . . .,h^'^''>^P\g{v{u))) over the real ze- 



(dcgps2,<T2-l), 
S2,o-2V-^' ■ ■ ■ •>PS2,(T2 

ros of p{u). 



4. Determine Sign going through the obtained list of realizable sign conditions. 

Proposition 10 Given geometric resolutions (psi,cri,t'5i,o-i,i, • • • ,i'Si,o-i,n) of Vsi,(7i ctnd 
{PS2,a2,vS2,(T2,i, ■■■, vs2,a2,n) ofVs2,a2, OLud the Thom encodings of a real root ofpsi,ai and Ps2,a2 
corresponding to points where the minimum of g on Vs-i,a^ Hi? and Vs2,a2 f^E is attained, Algo- 
rithm, ComparingMinimums compares these minimums within complexity 0{nD^^ ^^ + LDs^^s2 + 
-D^i,s2 log^ (-^51,^2)) where Z?si,s2 := max{Z:'sj , D^^ } with si = \Si\ and S2 = \S2\- 
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Proof. Step 1. This step is achieved within complexity 0(nZ)si,s2 log^(-Ds^_s2) loslos(-^si,s2)) 
as explained in Section [2j 

Steps 2, 3 and 4. Similar to Algorithm MinimumlnGeometricResolution Step 3(a), (b) and 
(c). The overall complexity of these steps is 0{nDl^ ^^ + LDs^^s2 + -^si S2 ^°S {^si,s2))- ^ 

We give now the main algorithm of the paper. 

Algorithm FindingMinimum 

INPUT: Polynomials /i, . . . , /„,, g € K[rEi, . . . , Xn] encoded by an sip of length L, an integer 
< I < 771 and an even integer d > deg(/j),deg(g). 

OUTPUT: A family {((pj,fj_i, . . . ,Uj^„),rj)}^^ where Z is a finite set and for every i ^ I, 
{Pi, Vi^i, . . . , Vi^n) is a geometric resolution in K[n] and tj G {—1, 0, l}™sPi jg ^j^e Thom encoding 
of a real root ^ of pi. 

1. Take a = (ai, . . . , a„) G K" at random and set ia '■= ctiXi + • • • + dnXn. 

2. S = {{S, a) \ S C {l,...,m} with < IS] < n and a € {+, -}^ with <7i = + for Z + 1 < 
i < m}. 

3. Take {Si,(7i) G S and remove it from S. 

4. {psi,ai,vsi,ai,i, ■ ■ ■ ,vsT_,ai,n) = GeometricResolution(/i(i G S'i),g, o-i,(i,£Q,). 

5. (Empty, r^^,^,,!, . . ■ ,Ts^,auk) = MiiiimuinIiiGeometricResolution(p5j,^-^,'(;5j,^-^,i, . . . ,vs^,cn,n, 
fi,- ■ ■ , fm,g,l,d). 

6. While Empty = True: 

(a) Discard (5i,(Ti), take a new (5i,(Ti) and remove it from S. 

(b) (P5i,ai,t^5i,f7i,i,---,'y5i,ai,n) = GeometricResolution(/i(i G S'ljjg, o-i,(i,^Q,). 

(c) (Empty, rs'^,,^-^,!, . . . .Ts^^a^^k) = MiniiiiumInGeometricResolution(p5'-^,<^-^, W5-^,^^,i, . . . , 
^S'l.o-i.n, fi,- ■ ■ , fm,g,l,d). 

8. While cS / 0: 

(a) Take (52, cr2) G S and remove it from S. 

(b) {PS2,o2i'"S2,a2,i^ ■ ■ ■ :Vs2,u2,n) = GeometricResolution(/i(i G S2),g,cr2,d,ea)- 

(c) (Empty, rsa, ,72,1, • • ■,rs2,a2,k) = MinimuinInGeometricResolution(p52,<^2i'"52,<T2,ii • • -,^82,^72,71, 
fi, ■ ■ ■ , fm,g,i,d). 

(d) If Empty = False : 

i. Sign = ComparingMinimums(p5^,^^,t;5^,^jj, . . ■ ,vsj,aun,PS2,a2,vs2,a2,i, ■ ■ ■ ,vs2,a2,n,g,d 

ii. If Sign = then 

2^ = 2^U{((pS'2,fT2)^S'2,CT2,l'--- )'^S'2,CT2,r!.)>TS'2,CT2,l)l ••• ' ((PS'2,CT2 > ^^2 ,f72,l ' ' ' ' ' ^S'2,<T2 ,n) , 7"52,(T2 ,fc') } " 
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iii. If Sign = 1 

^- ^= {((P52,CT2'^52,CT2,l'--- 'f^52,(T2,n),7"52,<T2,l)l ■■■, ((P52,(72 : '^52,<T2,1 ' • • • ^ '"S2,a2,n) , TS2,a2,k')} ' 

B. {Si,ai) = (82,(72). 

Now we prove Theorem 1. 

Proof of Theorem [IJ The correctness of the algorithm follows from the results in Sections [3] and 
m The complexity upper bound is obtained by adding the complexity bounds for the subroutines 
GeometricResolutionand MinimumlnGeometricResolution applied to every element in S, and 
ComparingMinimums applied successively to pairs of elements from S. D 
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